Effective Simulation of Quantum Entanglement Based on Classical Fields Modulated 

with Pseudorandom Phase Sequences 
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We demonstrate that n classical fields modulated with n different pseudorandom phase sequences 
can constitute a 2 n -dimensional Hilbert space that contains tensor product structure. By using 
classical fields modulated with pseudorandom phase sequences, we discuss effective simulation of Bell 
states and GHZ state, and apply both correlation analysis and von Neumann entropy to characterize 
the simulation. We obtain similar results with the cases in quantum mechanics and find that the 
conclusions can be easily generalized to n quantum particles. The research on simulation of quantum 
entanglement may be important, for it not only provides useful insights into fundamental features 
of quantum entanglement, but also yields new insights into quantum computation. 



The phenomenon of quantum entanglement is perhaps 
the most fascinating and important feature of quantum 
theory. On the one hand, quantum entanglement under- 
lies many of the most curious and controversial aspects 
of the quantum mechanical description of the world. On 
the other hand, quantum entanglement is widely appre- 
ciated as the essential ingredient of a quantum computer. 
[IHil . Quantum entanglement is first introduced by Ein- 
stein, Podolsky, and Rosen as most noticeable the EPR 
paradox [J] , which is at the origin of quantum nonlocal- 
ity. Bell proposed a remarkable inequality imposed by 
a local hidden variable theory, which enables an exper- 
imental test on the quantum nonlocality Some re- 
searchers have yet used classical approaches to simulate 
quantum entanglement theoretically and experimentally 
|6j-l8(. However, the simulation require an exponentially 
growth of physical resources with the number of quan- 
tum particles. This limitation could be traced back to 
a lack of a rigorous tensor-product structure in terms of 
subsystems js. [Toj. 

In this letter, we explore an effective simulation of 
quantum states based on classical fields modulated with 
pseudorandom phase sequences (PPSs). One essential 
difference between classical fields and quantum states is 
that n 2-state quantum particles could constitute a 2"- 
dimensional Hilbert space that contains tensor product 
structure, while n classical fields can not. In our pro- 
posal, by introducing the PPSs into the classical fields, we 
demonstrate that n classical fields modulated with n dif- 
ferent PPSs can also constitute a similar 2™-dimensional 
Hilbert space that contains a tensor product structure. 
In addition, we discuss classical simulation of Bell states 
and GHZ state, apply both the correlation analysis and 
von Neumann entropy to characterize them. We obtain 
similar results with the cases in quantum mechanics and 
find that these conclusions can be easily generalized to n 
quantum particles. 

The PPSs in our proposal derive from orthogonal pseu- 
dorandom sequences (PSs), which have been widely ap- 
plied to Code Division Multiple Access (CDMA) commu- 



nication technology as a way to distinguish different users 
111 Il2j . A set of PSs is generated by using a shift register 
guided by a Galois field GF(jp) that satisfies orthogonal, 
closure and balance properties [12j. In this letter, we con- 
sider an m-sequence of period N — 1 (N = p s ) generated 
by a primitive polynomial of degree s over GF(p) and 
apply it to 4-ary phase shift modulation, which has been 
a well-known modulation format in wireless and optical 
communications [ll| . A scheme is proposed to generate a 
PPS set S = {aW.AW,...^" 1 )} over GF(4). is 
an all-0 sequence and other sequences can be generated 
by using the method as follows: (1) given a primitive 
polynomial of degree s over GF(4) [13[, a base sequence 
of a length 4 s — 1 is generated by using Linear Feedback 
Shift Register; (2) other sequences are obtained by cyclic 
shifting of the base sequence; (3) by adding zeros to the 
sequences, the occurrence of any element equals to 4 S ~ 1 ; 
(4) mapping the elements of the sequences to [0,27r]: 
mapping 0, 1 mapping 7r/2, 2 mapping 7r, and 3 mapping 
3tt/2. 



We note the similarities between the Maxwell equation 
and the Schrodinger equation and the feasibility of using 
polarization or transverse modes of classical field to sim- 
ulate quantum states [14 16]. We first consider the two 
orthogonal modes (polarization or transverse) , which are 
denoted by |0) and |1) respectively, as the classical simu- 
lation of quantum bit (qubit) |0) and |1) @,[l6|. Thus any 
qubit state \ip) = a\0) + f3\l) can be simulated by a cor- 
responding classical mode superposition field expressed 
as \ip) = a|0)+j8|l), where |a| 2 + |/3| 2 = 1 (a,/3e C). 
Obviously, all the mode superposition fields could span a 
Hilbert space, where we can employ unitary transforma- 
tions to transform the mode status. 



Further, we discuss the classical fields modulated with 
multiple PPSs to simulate multiparticlc quantum system. 
Chosen two PPSs of A^ a ^ and from the set S, any 
two fields modulated with the PPSs can be expressed as 
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following 



= e lA(a> K|0) + /? a |l)). 
= e iA(b) K|0) + /3 b |l)). 



(1) 



According to the properties of the PPSs and the Hilbert 
space, we can define the inner product of the two fields 
and obtain the orthogonal property in our simulation, 
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1, a = b, 
0, a^b. 
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(a* a a b +/3* a j3 b ) 



(2) 



where A^ , A^ are the kth units of \( a \ A^, respec- 
tively. The orthogonal property supports to construct 
the tensor product structure of multiple fields. 

Assume two Hilbert spaces w and v spanned by the 
fields \4>a) and \ipb), any linear combinations of the el- 
ements in the direct product space of w ® v remain in 
the same space. We define the two orthogonal modes 
modulated with the PPS as the orthonormal bases 
of the space of w, expressed as |0 a ) = e lA<a) |0) and 
e lA< ' |1). Using the same notion, the orthonormal 
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\h) 

bases of the space of v are expressed as 1 0& ) = e lA '"' |0) 
and | lb) = e lA< ' |1). The four orthonormal bases are thus 
independent and distinguishable. Then the orthonormal 
bases of the direct product space of u>® v can be expressed 
as {|0„) ® |0 6 ) , |0.) ® |1 6 ) , |1.) ® \Q b ) , |1„) ® |l fe )}. Fur- 
therjWe can obtain the following tensor product proper- 
ties [l|: (1) for any scalar z, the elements \tp a ) , \tp b ) in the 
spaces of w and v, respectively, satisfy z (\ip a ) ® \ipb)) — 
(z \ip a )) ®\ipb) = \^a)®{z \ipb))] (2) in the space of w®v, 
the direct product of the combinations of elements equals 
to the combination of the direct products of elements, 

(\i>a) + Wa)) ® iW) + W b )) = \i> a ) ® W) + |V») ® + 

Wa) ® l^ft) + iV'a) ® iV'ft)- Using the same notion, we 
can construct a 2™-dimensional direct product space of 
\ip)® n = |^i) ® . . . ® \ip n ) by using the mode superposi- 
tion of n classical fields , . . . , \tp n ) modulated with n 
PPSs. 

Now we first discuss the correlation analysis of the clas- 
sical simulation of quantum entanglement and demon- 
strate the nonlocal correlation with Bell's inequality and 
equality criterion. As a contrast, we first study the prod- 
uct state of two classical fields expressed as Eq. ([1}, 
which corresponds to no entanglement in quantum me- 
chanics. In order to perform the correlation analysis, we 
need an orthogonal projection measurement. The result 
of projection measurement on \tp) can be obtained 
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FIG. 1: The scheme of the projection measurement. 

For convenience, the superposition coefficients a, ft are 
given by l/y/2, we obtain P (9) = cos (0). By using the 
projection measurement, we propose a correlation anal- 
ysis scheme in which the projection measurement of the 
basis |±) = |0) ± |1) is performed on each of the fields, as 
shown in Fig. [TJ P (9 a , k) = cos9 a and P (9 b , k) = cos 9b 
are the measurement results in the fcth sequence unit for 
the fields \ip a ) and \ip b ), respectively. The PPSs of A^ 
and A( fc ) do not affect the measurement results because 
they only contribute to the total phases. We can easily 
obtain the correlation function E(9 a ,9 b ) = cos 9 a cos 9b, 
which is similar to the case of quantum product state. 

Further, we consider that the modes |1) of the fields 
\ip a ) and \ipb) are exchanged by using a mode exchanger 
constituted by mode splitters and combiners, as shown 
in Fig. [5] We obtain the fields as following 



Wa) 
Wb) 



3 i7<a> |1' 



^ lb) IF 



(4) 



where the relative phase sequences (RPSs) 7™ = 
_ 7 (6) = A (b) _ A (a) j and 7 (a) + 7 (b) = 0. We obtain 

the results of the fields in the projection measurement 

P(9 a ,k) = cos (^+7^), P(0 b ,k) = cos (0 6 + 7^ 

where 7^,7^ are the fcth units of the RPSs ^ a \^ b \ 
respectively. Then we obtain the correlation function 

1 N _ _ 

E(9 a ,9 b ) = —Y, p (°a,k)P(9b,k) 



k=l 



COS ( 



(5) 



where C = 1/2 is the normalized coefficient. The corre- 
lation function is similar to the case of Bell state |\& + ). 
Therefore we consider the fields in Eq. (Q} as the classi- 
cal simulation of the Bell state |\E r+ ). We substitute the 
correlation functions above into Bell inequality (CHSH 
inequality) [I?} 



\B\ = \E(9 a ,9 b ) 
= 2^2 > 2, 



E(9 a ,9' b ) + E(9' a ,9' b )+E(9' a ,9 b )\ 

(6) 



where 9„ 



and 9' b are 7r/4, — 7r/4,0 and n/2, respec- 



(3) tively, when Bell's inequality is maximally violated. 
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FIG. 2: The mode exchange scheme for two classcial fields, 
where the gray blocks denote the mode splitters/combiners. 

Another Bell state l^ - ) differs from |\f r+ ) by tt phase. 
Similarly, we obtain the simulation of the Bell state 

I*") expressed as \ip' a ) = e lA<a) (jO) + e^ (a) |1)) /-y/2, 

W b ) = e 4A<b) (jO) + e ^ (b) +-) |1)) /y/2. By performing 
the transformation a x : |0) O |1) on \ip' b ) of the simula- 
tion of l^^), we obtain the simulation of the Bell state 

|$+) expressed as \tp' a ) = e lA<a> (jO) + e^ (a) |1)) /y/2, 
W b ) = e 4A<b) (e^ m |0) + |1)) /y/2, and of |$") ex- 
pressed as W a ) = e lA<a> (|0) + e^ (a) |1))/V2, W b ) = 

e ixW ^ e l ( 7<i>)+7r ) |0) + |1)) /y/2. Then their correlation 

functions Ey- (0 a ,0 b ) = - cos (0 a + b ), E^± {9 a ,9 b ) = 
±008(00. — Ob) are obtained. To substitute the correla- 
tion functions into Eq. (j6]), we also obtain the maxi- 
mal violation of Bell's inequality. The violation of Bell's 
criterion demonstrates the nonlocal correlation of the 
two classical fields in our simulation, which results from 
shared randomness of the PPSs. 

The nonlocality of the multipartite entangled GHZ 
states can in principle be manifest in a single measure- 
ment and need not be statistical as the violation of Bell 
inequality that relies on mean values (l8j . Preparing 
three classical fields \ip a ) , \ip b ) and \ip c ) similar to Eq. 
([1]), and circle exchanging the modes |1) of the fields, we 
obtain the fields as following 

Wa) = ^e* A(a) (|0) + e ^ (a) |l)), 
Wb) = ^"(lORe^), (7) 
Wo) = ^e iA(c) (|0) + ^ (c) |l)), 
where the RPSs 7 <» = - A( Q \ 7W = A( c ) - A^, 

7 (c) = \(a) _ A (c) j and 7 (a) + 7 (fc) + 7 (c) _ Q We ob _ 

tain the measurement results P (0 a , k) = cos (o a + J^^j , 

P {O b , k) = cos (e b + 7^) , P (9c k) = cos (0 C + 7 [. c) ) for 
the fields \ip' a ) , \ip' b ) and \ip' c ) in the projection measure- 
ment, respectively, and the correlation function 

1 N _ _ _ 

E(O a ,O b ,O c ) = —Y^P(Oa,k)P(O bl k)P(O c ,k) 

k=l 

= cos (O a + O b + Bo) , (8) 



where C — 1/4 is the normalized coefficient. If O a + 
O b + O c = 0, E(O a ,O b ,O c ) = 1. If O a + O b + O c = tt, 
E (0 a ,0b,0 c ) = — 1. This would manifest the nonlocal- 
ity in a single measurement for the simulation of GHZ 
states. It is noteworthy that the correlation function si- 
nusoidally oscillates with one of the phases O a , O b and O c , 
when the other phases fixed. However, the correlation 
function equals to zero if based on the correlation analy- 
sis of only two fields, because 7 ( - a - ) + + 7^ = and 
the sum of any two RPSs remains a pseudorandom se- 
quence. These results are completely similar to the case 
of GHZ states. 

Further, the simulation of GHZ state could be gener- 
alized to the case of n particles. Prepared n fields similar 
to Eq. (P) and circle exchanged the modes |1) of the 
fields, the RPSs satisfy 7M + . . . + 7^ = 0. We obtain 
the correlation function 

1 N _ _ 
E(9 1 ,...,9 n ) = Y, p {^)...P{6 n ,k) 

k=l 

= cos(0 1 + ... + 0„), (9) 

where P (0 t , k) = cos (Oi + 7^ is the result of the zth 
classical field at the kth sequence units in the projection 
measurement, and C — l/2 n ~ 1 is the normalized coeffi- 
cient. Using the same notion, the simulation of any other 
generalized GHZ states can be obtained and the corre- 
lation functions are also similar to the case of quantum 
mechanics. 

Now we continue to discuss our simulation in another 
view and try to apply von Neumann entropy as entangle- 
ment measure to characterize the classical simulation of 
Bell states and GHZ state. First we define a mode state 
in the direct product space of multiple classical fields as 
following 

»A< 3 > N 

m = e — E Wi) ® - ® K fc )) > a°) 

k=l 

where C is the normalized coefficient, A^' = ^"=1 ^ 
denotes the total phase sequence with the kth unit a£ , 
Wi°) denotes the mode superposition of the ith classical 
field at the kth sequence unit. Then, we introduce the 
density matrix formulation by using p — \fy) (^| and Eq. 
(fl"0|) . as following 

N 

p = E(^ S) Wl)®---®Wn)) 
fe=l 
N 

x E( elAia) (11) 
fe=i 

The entanglement of a partly-entangled pure state can 
be naturally parameterized by its von Neumann entropy 
of entanglement (l9j . Given a pure state p a b of two 
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subsystems a and b, we define the reduced density ma- 
trices p a = trb (pab) and p = tr a (p a b) for the states 
of a and b, where the partial trace has been taken 
over one subsystem, either a or b. Then the von Neu- 
mann entropy of the reduced density matrices is given by 
S = -tr (p a \og 2 p a ) = -tr (p b \og 2 pb). The quantity S 
ranges from zero for a product state to 1 for a maximally 
entangled pair of two-state particles, and < S < 1 for 
a partly entangled pair. 

We first discuss the entanglement measure of the 
product state of two classical fields expressed as 
Eq. ([T]). We can easily obtain the density matrix 
Pab = \ipa) (ipa\ ® l^b) (^i>|j and the von Neumann en- 
tropy S = -<r(p a log 2 (p a )) = ~tr(p b log 2 (pb)) = 
0. This indicates that no entanglement is involved 
between the two classical fields. We then discuss 
the simulation fields of Bell state |\& + ) expressed as 
Eq.Q. We could alter the expression as \^f) = 

e aW (|00) + |11)) /V2, where \q a q b ) , {q a ,b = 0, 1) are the 
orthonormal bases of the direct product space. It is 
worthwhile that the balance property J2k=i e%lh = 
of RPSs 7W results in many items of the direct 
product of |^) disappearing. Then we obtain p a b = 
(|00)(00| + |11) (11| + |00) (11| + |11) (00|)/2 and S = 
-tr (p a log 2 {pa)) = -tr (p&log 2 (p b )) = 1. It means the 
two fields of simulation of \^ + ) arc completely entangled. 
Using the same notion, we can obtain similar results for 
the simulation of other Bell states and GHZ state. Since 
the von Neumann entropy gives the same results for the 
simulation and quantum entangled states, we prove the 
validity of the simulation in a more rigorous way. 

It should be pointed out that the phase pseudorandom- 
ness provided by PPSs is different from the case of quan- 
tum mixed states. The latter results in decoherence of 
states and all coherent superposition items disappearing. 
Different from the decoherence, some coherent superposi- 
tion items remain in the mode state due to the constraints 
of the RPSs, such as 7W + 7^ = 0, 7^ + 7W + = 
for the simulation of Bell states and GHZ state, respec- 
tively. These remaining items makes the simulation of 
quantum entangled pure states possible. 

In this letter, we utilize the properties of PPSs to dis- 
tinguish classical fields that are overlapped in same space 
and time. A 2™-dimensional Hilbert space which contains 
tensor product structure is spanned by n classical fields 
modulated with PPSs. In our scheme, the resources re- 
quired are classical fields modulated with PPSs instead of 
optical/space modes. One classical field modulated with 
one PPS can simulate one quantum particle. It means 
that the amount of classical fields and PPSs grows lin- 
early with the number of quantum particles. According 
to the m-sequence theory, the number of PPSs in the set 
S equals to the length of sequences, which means that 
the time resource (the length of sequence) required also 
grows linearly with the number of the particles. Based on 



the analysis above, we conclude that one can efficiently 
simulate quantum entanglement with linearly growing re- 
sources by using our scheme. 

In summary, we have demonstrated a new scheme to 
simulate quantum entanglement and apply both the cor- 
relation analysis and von Neumann entropy to character- 
ize the simulation. We conclude that quantum entangle- 
ment can be efficiently simulated by using classical fields 
modulated with PPSs. The research on this simulation 
may be important, for it not only provides useful insights 
into fundamental features of quantum entanglement, but 
also yields new insights into quantum computation. 
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